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;_j ' We calculate the QCD moments of the pion structure function using Drell- 

Yan data on the quark distributions in the pion and a phenomenological model 
for the resonance region. The extracted higher twist corrections are found to 
\ be larger than those for the nucleon, contributing around 50% of the lowest 

moment at = 1 GeV^. 
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. Understanding the structure of the pion represents a fundamental challenge in QCD. As 

^ , the lightest qq bound state, the pion presents itself as somewhat of a dichotomy: on the 
Q \ one hand, its anomalously small mass suggests that it should be identified with the pseudo- 
' Goldstone mode of dynamical breaking of chiral symmetry in QCD, essential for describing 
^ the long-range structure and interactions of hadrons; on the other, high energy scattering 



experiments reveal a rich substructure which can be efficiently described in terms of current 
quarks and gluons. 
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Q \ The duality between quark and hadron degrees of freedom reveals itself in most spectac- 

' ular fashion through the phenomenon of Bloom-Gilman duality in inclusive deep-inelastic 
scattering. Here the inclusive F2 structure function measured at low hadron final state mass 
^ . W , in the region dominated by low-lying resonances, has been found in the case of the proton 
[1,2] to follow a global scaling curve which describes high W data, to which the resonance 
structure function averages. Furthermore, the equivalence of the averaged resonance and 
scaling structure functions for each prominent resonance region separately suggests that the 
resonance-scaling duality also exists to some extent locally. 

Within QCD, the appearance of Bloom-Gilman duality for the moments of structure 
functions can be related through the operator product expansion to the size of high twist 
corrections to the scaling structure function [3] . The apparent early onset of duality for the 
proton structure function indicates the dominance of single-quark (leading twist) scattering 
even at low momentum transfers. It is not a priori clear, however, whether this is due to 
an overall suppression of coherent effects in inclusive scattering, or because of fortuitous 
cancellations of possibly large corrections. Indeed, there are some indications from models 
of QCD that the workings of duality may be rather different in the neutron than in the 
proton [4,5], or for spin- independent and spin-dependent structure functions. Given that 
Bloom-Gilman duality is empirically established only for baryons (specifically, the proton), 
while the application of theoretical models is generally more straightforward in the meson 
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sector, a natural question to consider is whether, and to what extent, duahty manifests itself 
phenomenologically for the simplest qq system in QCD — the pion. 

In this note we report the first analysis of the role of resonances in the QCD moments of 
the pion structure function, and obtain a first estimate of the size of higher twist corrections 
to the scaling contribution. Similar analyses for the nucleon have been made in Refs. [6-8] . 
According to the operator product expansion in QCD, at large the moments of the pion 
F2 structure function, 

M„(Q2)= f'dxx^-' F^{x,Q') , (1) 
Jo 

can be expanded as a power series in with coefficients given by matrix elements of 

local operators of a given twist r, 

00 / 1 \ 

Mn{Q')^T.AUMQ'))[^Q^j . (2) 

Here the leading twist r = 2 term Aq corresponds to free quark scattering, and is responsible 
for the scaling of the structure functions (modulo perturbative as{Q^) corrections). The 
higher twist terms ^">2 represent matrix elements of operators involving both quark and 
gluon fields, and are suppressed by powers of 1/Q'^. The higher twist contributions refiect 
the strength of nonperturbative QCD effects, such as multi-parton correlations, associated 
with confinement. 

Note that the definition of M„((5^) includes the elastic contribution a.t W = m^, corre- 
sponding to a; = Q'^/{W'^ — rri^ + Q^) = 1, which is given by the square of the elastic pion 
form factor, Ft^{Q'^), 

f;^'\x = 1, Q') = 2m^u {F^{Q')f S{W' - ml) . (3) 

Although negligible at high Q^, the elastic contribution has been found to be important 
numerically at intermediate for moments of the nucleon structure function [6-10]. In 
Eq. (1) we use the Cornwall- Norton moments rather than the Nachtmann moments, which 

are expressed in terms of the Nachtmann scaling variable, = 2x/{l + \Jl + 4x^m^/(5^), 
and include effects of the target mass. Because of the small value of m^, the difference 
between the variables x and ^, and therefore between the x- and ^-moments, is negligible 
for the pion. 

The pioneering analysis of De Rujula et al. [3] (see also Ref. [6]) showed that the size of 
the higher twist matrix elements directly governs the onset of quark-hadron duality. Namely, 
there is a region of n and in which the moments of the structure function are dominated 
by low mass resonances, where the higher twist contributions are neither overwhelming 
nor negligible. For example, even though there are large contributions from the resonance 
region (PF < 2 GeV) to the n — 2 moment of the proton F2 structure function (~ 70% at 
Q"^ = 1 GeV^), the higher twist effects are only of the order 10-20% at the same [6]. The 
question we address here is whether there is an analogous region for the pion, where the 
resonance contributions are important, but higher twist effects are small enough for duality 
to be approximately valid. 
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Of course, strictly speaking the distinction between the resonance region and the deep 
inelastic continuum (DIS) is somewhat arbitrary, as can be illustrated in the large A^^ limit of 
QCD. There the final state in deep inelastic scattering from the pion is populated by infinitely 
narrow resonances, even in the Bjorken limit, while the structure function calculated at the 
parton level produces a smooth, scaling function [4,11]. Phenomenologically the spectrum 
of the excited states of the pion is expected to be rather smooth sufficiently above the p 
mass, for ^ 1 GcV. Contributions from the excitation of heavier mesons are not expected 
to be easily discernible from the DIS continuum — the ai meson, for instance, appears at a 
mass W ~ 1.3 GeV, and has a rather broad width (~ 350-500 MeV). 

The pion structure function has been measured in the t^N Drell-Yan process [12-14] 
over a large range ofx, 0.2^a;^l, and for typically ^ 20 GeV^. It has also been 
studied at HERA in semi-inclusive DIS at very low x and high W [15]. No data exist on 
F2 at low W , however, in the region where mesonic resonances would dominate the cross 
section. The spectrum could in principle be reconstructed by observing low t neutrons at 
low W produced in the semi-inclusive charge-exchange reaction, ep — > enX, where t is the 
momentum transfer squared between the proton and neutron. In the absence of such data, 
to obtain a quantitative estimate of the importance of the resonance region, we model the 
pion spectrum at low W in terms of the elastic and p pole contributions, on top of the 
DIS continuum which is evaluated by evolving the leading twist structure function to lower 
Q"^. The latter can be reconstructed from parameterizations [16-18] of leading twist quark 
distributions in the pion obtained from global analyses of the pion Drell-Yan data. In this 
work we use the low fit from Ref. [16], which gives the leading twist parton distributions 
in the pion for > 0.25 GeV^, although our conclusions do not change with the use 
of other parameterizations [17,18]. For the elastic contribution we use a parameterization 
of the world's data [19] which interpolates smoothly between the perturbative QCD and 
photopro duct ion limits [20]. 

The p contribution is described by the vr — > p transition form factor, F^^plQ'^), and is 
expected to fall as l/Q*^ at large (compared with for Ft^^Q"^)). Since there is no 

empirical information on Fj^plQ"^), we consider several models in the literature, based on a 
relativistic Bethe-Salpeter vertex function [21], a covariant Dyson- Schwinger approach [22], 
and light-cone QCD sum rules [23]. These represent a range of ~ 100% in the magnitude of 
FnpiQ'^) over the region of covered in this analysis. The calculation of Ref. [23] gives a 
somewhat smaller result than in Rcfs. [21,22], which give a similar magnitude for F-^p. The 
difference between these can be viewed as an estimate of the uncertainty in this contribution. 

In Fig. 1 we plot the contributions to the moments of the pion structure function from 
the "resonance region", M^^{Cp), as a ratio to the total moment, for n = 2, • • • , 10. The 
resonance region here is defined by < W^rcs = 1 GcV, corresponding to restricting the 
integral in Eq. (1) to the range x^-^s = Q"^ /(Wrcs ~ + Q^) ^'^ ^^c elastic point at x = 1. 
For the n = 2 moment the low W region contributes as much as 50% at = 2 GeV^, 
decreasing to < 1% for > 10 GeV^. Higher moments, which are more sensitive to the 
large x region, subsequently receive larger contributions from low W. The n = 10 moment, 
for example, is virtually saturated by the resonance region at = 2 GeV^, and still has 
some 40% of its strength coming from W < 1 GcV even at = 10 GeV^. 

The large size of the resonance contributions suggests that, at a given scale Q^, higher 
twist effects play a more important role in the moments of the pion structure function than 
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FIG. 1. Contributions to moments of the pion structure function from the resonance region, 
W < Wres = 1 GeV, relative to the total. 

in the case of the nucleoli. The lowest (n = 2) moment of F2 is displayed in Fig. 2, including 
the leading twist and elastic contributions to M2{Q^). The leading twist component, 



is expressed (at leading order in as{Q^)) in terms of the twist-2 quark distributions in the 
pion, q'^{x,Q'^). The leading twist contribution is dominant at > 5 GeV^, while the 
deviation of the total moment from the leading twist at lower indicates the increasingly 
important role played by higher twists there. In particular, while negligible beyond ^ 
4 GeV^, the elastic contribution is as large as the leading twist already at ~ 1 GcV^. The 
contribution from the tt — > p transition is more uncertain, and the band in Fig. 2 represents 
the total moment calculated using different models [21-23] of Fj^p{Q'^). However, while the 
current uncertainty in this contribution is conservatively taken to be ~ 100%, doubling 
this would lead to a modest increase of the band in Fig. 2. Uncertainty arising from poor 
knowledge of the leading twist distributions at small x [16-18] is not expected to be large. 

The higher twist part of the moments can be defined as the difference between the total 
moment M^iQ"^) and the leading twist contribution in Eq. (4), which includes a term arising 
from target mass corrections. 



Although nonperturbative effects can in principle mix higher twist with higher order effects 
in q;^, rendering the formal separation of the two ambiguous [24] (the perturbative expansion 




(4) 



(5) 



4 



0.6 




\ 



leading twist 



\ 











elastic 

1 """^ — < 

2 4 



6 



J. 



8 



10 



(2^(GeV2) 



FIG. 2. Lowest (n = 2) moment of the pion structure function. The leading twist (sohd) and 
elastic (dashed) contributions are shown, and the shaded region represents the total moment using 
different models for the tt ^ p transition. 

itself may not even be convergent), by restricting ourselves to the region of in which the 
l/Q^ term is significantly larger than the next order correction in a^, the ambiguity in 
defining the higher twist terms can be neglected [6]. Because the target mass correction, 
M™((5^), which is formally of leading twist, is proportional to m^/Q^, its contribution 
will only be felt when ~ m^, which is far from the region where the twist expansion is 
expected to be valid. The higher twist contribution to the n = 2 moment is plotted in Fig. 3 
as a ratio to the total moment, as a function of Q^. The band represents an estimate of the 
uncertainty in the tt — > p transition form factor, as in Fig. 2. The higher twist contribution 
is as large as the leading twist at = 1 GeV^, is around 1/3 at = 2 GeV^, and vanishes 
rapidly for > 5 GeV^. 

The higher twist contribution at ~ 1 GeV^ appears larger than that found in sim- 
ilar analyses of the proton F2 [6] and gi [8] structure functions. This can be qualitatively 
understood in terms of the intrinsic transverse momentum of quarks in the hadron, (k^), 
which typically sets the scale of the higher twist effects [3,6,25,26]. By analyzing the x — > 1 
dependence of the measured n'^fi~ pairs produced in vriV collisions, and the angular distri- 
bution at large x, the E615 Collaboration [14] indeed finds the value (k'^) = 0.8 ± 0.3 GeV^ 
within the higher twist model of Ref. [25], which is larger than the typical quark transverse 
momentum in the nucleon (C(500 MeV)). The implication is that duality would therefore 
be expected to set in at larger for the pion than for the nucleon. 

In summary, we have evaluated moments of the pion structure function, and studied in 
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FIG. 3. Higher twist contribution to the n = 2 moment of the pion structure function, as a 
ratio to the total moment. The band indicates the uncertainty due to the model dependence of 
the TT — > p transition form factor. 

particular the role of the resonance region. Making the reasonable assumption that the low 
W resonant spectrum is dominated by the elastic and tt — > p transitions, we have presented 
a first quantitative estimate of the size of higher twist contributions. For the lowest moment 
of F2 we find that the resonance region (W < 1 GeV) contributes ~ 50% of the strength 
at ~ 2 GeV^, dropping to below 10% for ^ 5 GeV^. The elastic component, while 
insignificant for ^ 3 GeV^, is as large as the leading twist contribution at ~ 1 GeV^. 
The higher twist corrections to the n = 2 moment amount to ~ 50% at = 1 GeV^, 
~ 30% at = 2 GeV^, and become neghgible beyond ~ 6 GeV^. 

Uncertainties in these estimates are mainly due to the poor knowledge of the inclusive 
pion spectrum at low W, which limits the extent to which duality in the pion can be tested 
quantitatively. Only the elastic form factor has been accurately measured to ^ 2 GeV^, 
although at larger it is poorly constrained. The inclusive pion spectrum can be extracted 
from data from the semi-inclusive charge- exchange reaction, ep — > enX, at low t, for instance 
with CLAS at Jefferson Lab. This could also allow one to determine the individual exclusive 
channels at low W. In addition, a Rosenbluth separation would allow the transverse and 
longitudinal structure functions of the pion to be extracted. 
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